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Abstract 

In  phase  and  frequency  measurements,  the  measured  phase 
and  frequency  are  not  the  true  phase  and  frequency  but  the  ones 
which  are  disturbed  by  noises,  due  to  the  effects  of  the  noise 
processes.  In  this  paper,  we  discussed  the  effects  of  three  noise 
processes,  i.e.,  White  PM,  White  FM  and  Random  Walk  FM,  on  the 
estimations  of  phase  and  frequency.  It  is  indicated  that  the 
properties  of  these  two  estimations  are  very  different.  In  phase 
estimation,  the  error  can  be  reduced  by  properly  selecting  suitable 
smoothing  length  NT  and  smaple  interval  T,  But  in  frequency 
estimation,  the  error  cannot  be  reduced  arbitarily  by  means  of 
improving  estimator  or  measurement  equipment.  The  precison  of 
frequency  is  limited  by  the  intrinsic  noises  in  the  clock. 


I.  Introduction 

In  phase  and  frequency  measurements,  the  measured  phase  and  frequency 
are  not  the  true  phase  and  frequency  but  the  disturbed  ones,  due  to  the  effects 
of  the  noise  processes.  Before  an  atomic  time  scale  being  computed,  the  phase 
and  frequency  of  each  clock  should  be  measured.  So  it  is  necessary  to  give 
proper  estimations  of  the  true  phase  and  frequency. 

In  this  paper,  we  discussed  the  effects  of  three  noise  processes,  i.e.,  White 
PM,  White  FM  and  Random  Walk  FM,  on  the  estimations  of  the  phase  and 
frequency.  It  is  indicated  that  the  properties  of  these  two  estimations  are  very 
different. 

II.  The  Noise  Model 

The  basic  model,  so  called  Power-law  Model,  is  presented  in  frequency 
domain[l],  i 

Sy(£)  ~h_af'*+h_1f-1+h0+h1f+hafa  ( 1 ) 

In  time  domain,  the  Dynamic  Model  is  more  convenient[2], 

m)=xU)+r{ic) 

x{k+l)  ~x{k)  +y(Jc)  T+q{k)  (2) 

y(k+l)-y(k)+f{k) 

where  r(k)  are  phase  modulation  noise  processes  which  contain  White  PM  and 
Flicker  PM;  q(k)  are  frequency  modulation  noise  processes  including  White  FM, 
Flicker  FM  and  Random  Walk  FM;  and  f{k)  is  Random  Walk  FM.  The  sampling 
interval  is  T. 

The  two  flickers  are  not  considered  in  the  following  disscusion,  because  the 
other  three  noise  processes  in  different  levels  are  enough  to  demenstrate  the 
properties  of  thfe  phase  and  frequenbcy  estimations.  So  r(k)  represents  the  White 
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including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

DEC  1992 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-1992  to  00-00-1992 


4.  TITLE  AND  SUBTITLE 

The  Limit  of  Frequency  Estimation 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Shaanxi  Astronomical  Observatory, Academia  Sinica,P.  0.  Box  18 
Lintong, 710600  Shaanxi,  P.R.  China, 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA267301.  24th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Applications  and  Planning 
Meeting,  McLean,  VA,  1-3  Dec  1992 


14.  ABSTRACT 

see  report 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

5 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


PM  only, 


Dr(k)2=Er(k )2*  -  1  -h1fh  (3) 

(2«)2 

and  q(k)  represents  the  White  FM  and  the  Random  Walk  FM, 

Dq(k)  2=Eq{k)2  =  ^-h0T^~  (2  n)2h.2T3  (4) 

2  6 

III.  Phase  Estimation  (finite- memory  smoothing) 

For  a  set  of  2AN-1  phase  measurements  of  T{k+i)(i=-Ni...10i...,N)i  the  linear 
estimation  of  Mi)  is 

s 

£(k)  =  J2  CiTik+i)  (5) 

There  is  a  set  of  cofficients,  Ci(i=-JVt...10,...,./V),  which  enables  the  estimation 
error  to  be  minum, 

£[m)-J?(Jt)3^  min  (6) 

But  we  do  not  care  about  the  optimal  solution  here,  because  the  properties  of  the 
estimation  are  more  important.  For  convenience,  a  simplified  smoother,  average 
smoother,  ip  used  here.  That  is, 

Ci=  Zifel  i-N.-'O'-'ir  (7) 

and 

*(k)  - 


If  T  is  not  large,  the  effect  of  Random  Walk  FM,  f[k)t  can  be  ignored  in  phase 
estimation.  It  is  easy  to  show  that 

i 

x(k+i)  =x(k) q(k+l-l)  (9) 

So  the  estimated  phase  can  be  expressed  as  true  phase  plus  a  serier  of  noise 
processes. 


x 


2K+1 


£  m+i) 


2WA 


V  x{k+i)  + 


2N+1 


K 

£ 


(8) 


r{k+i) 


(10) 


- qik~d) 

t"2&r£i 

According  the  assumption  of  independence  between  noise  processes,  the 
estimation  error  is 
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E[x{k)  -*<*}] 2 

-zhr'-T^: -*A+  rm£i)b°T 


£  [N-i+l]  imEq2  (k) 


(11) 


2N+1  (2ir) 2 

»JL- — h. 


'T"* 

This  is  the  property  of  phase  estimation  that  a  small  NT,  smoothing  length, 
can  reduce  the  effect  of  White  FM  and  also  when  NT  is  fixed,  increasing  N,  i.e., 
decreasing  T,  can  reduce  the  effect  of  White  PM.  Therefore,  a  suitable  selection 
of  N  and  T  will  make  the  error  of  phase  estimation  small  enough. 

IV.  Frequency  Estimation- 1  (Finite-memory  smoothing) 

We  assume  the  measured  phase  has  been  carefully  filtered  before  the 
frequency  estimation  being  done.  Thus  the  phase  noise  can  be  ignored  and  the 
noise  model  of  Eq.(l)  is  rewritten  as 


where 


d(k)  my(k)  +-iqr(Jc) 
y{k+l)  ~y{k)  +f{k ) 


d(*U-i;[x(it+l)-JC<Jc)] 


(12) 


(13) 


For  a  set  of  2N+1  frequency  measurements  of  cf(/c+i)(i=-~Ar,...,0,...»jV},  the 
linear  estimation  of  y(t)  is 


*<*)-£  Ctdik+i) 


(14) 


We  also  consider  an  average  smoother,  similar  to  the  case  in  phase 
estimation.  That  is, 


'w"ssr£<«**" 

■ ski fl 1W)  w-1*” 

tRr£«*i-1>  W-U1) 


(15) 


The  estimation  error  is 
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E[y[k) -f[k))z 

1  o  N 

■aKT  [T?h«+l  <2,t,^I1+i^r-i  (2'>lh-‘T 


2T 

•~NTh0+it1 2h_2NT 


(16) 


It  can  be  seen  from  Eq.(H)  and  Eq.(17)  that  estimations  of  phase  and 
frequency  ate  very  different.  In  phase  estimation,  the  two  noise  effects  can  be 
reduced  simultaneously  by  small  NT  and  large  N.  But  in  frequency  estimation 
small  NT  may  reduce  the  effect  of  h_2  but  the  effect  of  h0  -will  increase,  while 
large  NT  will  lead  high  h_2  and  low  h~.  Therefore,  the  precision  of  frequency 
estimation  cannot  be  improved  arbitarily  by  change  N  or  T  for  certain  noise 
levels. 

According  to  Eq.(17),  we  may  choose  a  suitable  NT  which  can  make  the 
estimation  error  as  small  as  possible, 


[NT) 


min' 


2  It  \ 


U7) 


Correspondent  error  is 


E[y{k)  -5?Uc)  ]  siin** 


(18) 


V.  Frequency  Estimation-II  (Kalman  Filtering) 


For  the  noise  model  shown  in  Eq,(12),  Kalman  filter  may  be  considered  as 
an  optimal  estimator.  Obviously,  the  system  of  Eq.(12)  fulfills  the  conditions  of 
observability  and  controllability,  so  the  filter  is  stable.  When  the  system  is  in 
steady  state,  the  frequency  estimation  is 

?(k+l/k+l)  =j?(Jc/Jc)  +Kld(k+1)  -?(k/k)  ]  < 19  > 

and  the  estimation  error  is 

Eiy[k)-?[k/k)]2=P+  (2°) 

In  steady  state,  the  filtering  error  P*,  the  prediction  error  P  and  the  gain  K 
fulfill  the  Riccati  Equation  as  below, 


P'=P++F 

IC«p-<p-*-i-0)  -1 
T2 

P*=(1-F)P~ 

where  F=Dl*(k)  and  Q=D(q(k)/T)2,  P  is  solved  as 


(21) 


1  (-F+ 

2 


Fa+4F-^-  } 

j>2 


2h_2T+ 

There  is  a  Tmin  which  makes  P  be  minium, 


(22) 


n2hzh_2+~n*hl2T2 
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(23) 


T  ml. 
mln 

The  correspondent  estimation  error  is 


9  An 


28  h_ 


B[y(k)  -?(k/k)] 


(24) 


Because  Kalman  filter  is  an  optimal  estimator,  any  other  estimator  will  not 
be  better  than  it.  Therefore,  Eq.(25)  can  be  regarded  as  the  limit  of  the 
frequency  estimaton.  This  limit  is  related  to 'the  intrinsic  noises  in  the  clock,  not 
the  measurement  equipment.  In  general  speaking,  the  precision  of  frequency 
estimation  is  limited  by  the  noise  levels  in  the  clock,  which  cannot  be  improved 
arbitarily  by  improving  measurement. 


VI.  Conclussion 


We  have  analyzed  the  effects  of  clock  noises  on  the  phase  and  frequency 
estimations,  respectively.  In  phase  estimation,  the  error  can  be  reduced  by 
properly  selecting  suitable  smooth  length  NT  and  smaple  interval  T.  But  in 
frequency  estimation,  the  error  cannot  be  reduced  arbitarily  by  means  of 
improving  estimator  or  measurement  equipment.  The  precison  of  frequency  is 
limited  by  the  intrinsic  noises  in  the  clock. 
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